Search for: All records

The most general bound binary black hole system has an eccentric orbit and precessing spins. The detection of such a system with significant eccentricity close to the merger would be a clear signature of dynamical formation. In order to study such systems, it is important to be able to evolve their spins and eccentricity from the larger separations at which the binary formed to the smaller separations at which it is detected, or vice versa. Knowledge of the precessional evolution of the binary’s orbital angular momentum can also be used to twist up aligned-spin eccentric waveform models to create a spin-precessing eccentric waveform model. In this paper, we present a new publicly available code to evolve eccentric, precessing binary black holes using orbit-averaged post-Newtonian (PN) equations from the literature. The spin-precession dynamics is 2PN accurate, i.e., with the leading spin-orbit and spin-spin corrections. The evolution of orbital parameters (orbital frequency, eccentricity, and periastron precession), which follow the quasi-Keplerian parametrization, is 3PN accurate in the point particle terms and includes the leading order spin-orbit and spin-spin effects. All the spin-spin terms include the quadrupole-monopole interaction. The eccentricity enhancement functions in the fluxes use the high-accuracy hyperasymptotic expansions from Loutrel and Yunes []. We discuss various features of the code and study the evolution of the orbital and spin-precession parameters of eccentric, precessing binary black holes. In particular, we study the dependence of the spin morphologies on eccentricity, where we find that the transition point from one spin morphology to another can depend nonmonotonically on eccentricity, and the fraction of binaries in a given morphology at a given point in the evolution of a population depends on the instantaneous eccentricity. 

Subjected to the tidal field of its companion, each component of a coalescing binary suffers a slow change in its mass (tidal heating) and spin (tidal torquing) during the inspiral and merger. For black holes, these changes are associated with their absorption of energy and angular momentum fluxes. This effect modifies the inspiral rate of the binary, and consequently, the phase and amplitude of its gravitational waveform. Numerical relativity (NR) waveforms contain these effects inherently, whereas analytical approximants for the early inspiral phase have to include them manually in the energy balance equation. In this work, we construct IMRPhenomD_Horizon, a frequency-domain gravitational waveform model that incorporates the effects of tidal heating of black holes. This is achieved by recalibrating the inspiral phase of the waveform model IMRPhenomD to incorporate the phase corrections for tidal heating. We also include corrections to the amplitude, but add them directly to the inspiral amplitude model of IMRPhenomD. First we demonstrate that the inclusion of the corrections, especially in the phase, confers an overall improvement in the phase agreement between the analytical inspiral model (uncalibrated SEOBNRv2) and NR data. The model presented here is faithful, with less than 1% mismatches against a set of hybrid waveforms (except for one outlier that barely breaches this limit). The recalibrated model shows mismatches of up to ∼14% with IMRPhenomD for high mass ratios and spins. Amplitude corrections become less significant for higher mass ratios, whereas the phase corrections leave more impact—suggesting that the former is practically irrelevant for gravitational wave data analysis in Advanced LIGO (aLIGO), Virgo and KAGRA. Comparing with a set of 219 numerical relativity waveforms, we find that the median of mismatches decreases by ∼4% in aLIGO zero-detuned high power noise curve, and by ∼1.5% with a flat noise curve. This implies a modest but notable improvement in waveform accuracy.